Fermionic atoms trapped in one-dimensional optical superlattice with harmonic 

confinement 
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We study the ground-state properties of spin- 1/2 fermionic atoms confined in a one-dimensional 
optical superlattice with harmonic confinement by using the density matrix renormalization group 
method. For this purpose, we consider an ionic Hubbard model that has superlattice potentials 
with 2-site periodicity. We find that several different types of insulating regimes coexist even if the 
number of atoms at each site is not an integer, but its average within the unit cell is an integer 
or half integer. This is contrasted to the coexisting phase of the metallic and Mott-insulating 
regimes known for the ordinary Hubbard model in an optical lattice. The phase characteristics are 
elucidated by investigating the profiles of the atom density, the local density/spin fluctuations, the 
double occupation probability and the spin correlations in detail. 

PACS numbers: 03.75.Ss, 05.30.Fk, 34.50.-s, 71.30.+h 



I. INTRODUCTION 

Experimental techniques for manipulating ultracold 
atoms have made great progress since the successful re- 
alization of atomic Bose-Einstein condensates (BECs) 
EJ. Optical lattices, formed by a standing wave of laser 
light, are providing the ideal stages for an experimen- 
tal investigation of the fundamental many-body prob- 
lems in condensed matter physics via ultracold atomic 
gases [S IS [J 13 A recent series of experimental studies 
on bosonic Mott-insulators [S , IS El ES El El have 
clearly demonstrated this feature, namely the realization 
of quantum simulators. 

A noteworthy advantage of atomic gases over con- 
densed matter such as solids or liquids is that their exper- 
imental parameters are highly controllable. Both density 
and temperature of atomic gases are fully controlled in 
the process of evaporative cooling |l| . The depth of the 
optical lattice potential is controlled by the light inten- 
sity, which leads to the precise alternation of atom tun- 
neling between adjacent sites 0,0- Furthermore, we can 
create one-, two-, and three-dimensional optical lattices 
depending on the laser beam configuration with the use 
of interference effects 0, Even the interactions be- 
tween the atoms can be widely tuned with the help of 
Feshbach resonances El- Therefore, in the experiments, 
it is possible to precisely change the important system 
parameters such as density, temperature, tunneling, di- 
mensionality of lattices, interaction strength and so forth. 
Ultracold atoms in optical lattices allow us to investigate 
the quantum many-body effects in a well-controlled man- 
ner. 

Experiments with quantum degenerate fermionic gases 
have already clarified the intri guing ph ysics of such phe- 
nomena as a molecular BEC [Tj,[ia,ll3 an d the crossover 
between a BEC and a Bardeen-Cooper-Schrieffer (BCS) 



superfluid [TtI ITsl filj. l2p| |. Ultracold fermionic atoms 
trapped in optical lattices are now stimulating both ex- 
perimental and theoretical interest. This system is well 
captured by the Hubbard model Hamiltonian |2l| , which 
is widely discussed with respect to strongly correlated 
electron systems. Several theoretical analyses have pre- 
dicted the possibility of the Mott metal-insulator tran- 
sition (MMIT) in this system [U H H |H 
Althou gh e xperimental investigations are still under way 
[271 1281 129L l30| , the recent demonstration of a fermionic 
band insulator in three-dimensional optical lattices |29| 
has convinced us that MMIT will be observed in near 
future. A precise comparison of the theory and experi- 
ments will provide us with a deeper understanding of the 
quantum phase transition in fermionic systems. 

The superlattice geometry of optical lattices, i.e., op- 
tical superlattice, was successfully demonstrated as an 
advanced manipulation technique in a recent report [3]| . 
This technique has been introduced to control the dis- 
tribution of atoms on lattices more functionally and will 
play a key role in the implementation of quantum com- 
puting with ultracold atoms 0, |^| . On the other hand, 
in condensed matter physics, it is known that lattice 
distortions generate a variety of many-body effects in 
strongly correlated electron systems such as the charge- 
transfer or gan ic materials and ferroelectric perovskites 
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The ultracold fermionic atoms in optical superlattices 
01 can he regarded as a quantum simulator for these 
materials. Our present study is highly motivated by these 
important potentialities of optical superlattices. 

In this paper, we investigate the ground state prop- 
erties of S=l/2 fermions trapped in a one-dimensional 
optical superlattice with 2-site periodicity. The external 
harmonic confinement inherent in the real experiments 
|S| is further considered in our model. We show that a 
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FIG. 1: Potential profiles of (a) 2-site periodic superlattice 
and (b) with a confining potential. Vsl is the difference be- 
tween the potential energies of odd and even sites. 



new coexistence phase appears consisting of three distinct 
insulating regions caused by the superlattice potentials, 
which is in contrast to the ordinary Hubbard model in 
the optical lattice. This is elucidated via a detailed anal- 
ysis of the local density of atoms and spin correlations by 
using the density matrix renormalization group method. 
Furthermore, we find unique spin correlations between 
the insulating regions when these regions are separated 
by a metallic region. 

This paper is organized as follows. In the next section 
we briefly mention the model and the method. We then 
show the density profiles of fermions and several den- 
sity/spin correlation functions to clarify the zero tem- 
perature properties. The final section provides a brief 
summary. 

II. MODEL AND METHOD 

We consider a gas of fermionic atoms embedded in an 
optical superlattice with an alternating superlattice po- 
tential and a confining parabolic potential (see Fig. 1(b)), 
which is described by the Hubbard Hamiltonian, 

H = -t ^](ct ff Cf + i ;ff + h.c.) + U n^n %i 

i,u i 

+ ]T^7W+I4 5>-iV s /2) 2 71 m (1) 
z,cr ia 

where ct , ci a and ni a — c\ a Ci,j are respectively the 
fcrmion creation, annihilation and number operators rel- 
evant to the site labeled i with spin cr(=t, J,). Here, U 
is the on-site repulsion (U > 0), V c the curvature of the 
parabolic confining potential, Vi the local potential at site 
i and t the hopping amplitude between nearest neighbor 
sites. 

We define the local potential Vi with 2-site periodicity 

as 

K = tfcL{mod(i+l,2)} (2) 
where Vs l is the difference between the potential energies 



of odd and even sites (see Fig. 1(a)). We denote the total 
number of fermions (lattice sites) as Nf (N s ). 

There are several kinds of numerical methods that can 
be employed with our fermionic superlattice Hubbard 
model. The numerical diagonalization method can be 
applied straightforwardly to fermionic models, but it re- 
stricts our analysis to rather small systems with a typ- 
ical lattice size of 20 sites, which is not large enough 
for our purpose. The quantum Monte Carlo simulations 
method is an efficient way to deal with quantum systems 
at finite temperatures, and this approach has been ap- 
plied successfully to bosonic Hubbard models to elucidate 
the transition between the superfluid and Mott insulat- 
ing phases in optical lattices 0i Oil EH ■ This method, 
however, suffers from negative-sign problems in fermionic 
cases, which makes it difficult to obtain reliable results 
at low temperatures. The density matrix renormaliza- 
tion group (DMRG) |5lL IH^ is another powerful numeri- 
cal method with which we investigate quantum systems 
such as fermionic models with extremely high accuracy, 
although it is restricted to one-dimensional cases. This 
technique allows us to deal with a very large system as 
required for the analysis of optical lattices. In this pa- 
per we therefore exploit a variant of the DMRG (the so- 
called finite-system DMRG) 51] to study our fermionic 
Hubbard model in an optical superlattice at zero temper- 
ature. The essence of the idea is that we first diagonalize 
a given small system, and then enlarge the system size 
step by step via a renormalization procedure by retaining 
relevant low-energy states which have large eigenvalues 
of the density matrix. This simple procedure is known 
to give precise results for one-dimensional quantum sys- 
tems. This method has already been reported in detail 
|5lLl52| . In our nonuniform model with harmonic confine- 
ment, the numerical calculation can be done in parallel to 
the ordinary case without confinement. The only point 
we have to pay attention to is that the superlattice poten- 
tial should be properly included in each renormalization 
step. 

We restrict our discussions to the singlet ground state 
with an even number of fermions, which are tractable by 
the DMRG method [5l|,E2- As far as the local density 
and its variance are concerned (see Sec. Ill A and B), the 
even-number case is sufficient to discuss the characteristic 
properties. A relevant difference between the even and 
odd cases may appear in the spin correlations (Sec. Ill 
C) since a free S=l/2 spin remains unscreened in the 
system for the odd case. The free spin is not important 
physically since it does not affect the global nature of the 
system, but rather acts as a magnetic impurity. 



III. RESULTS 

A. Density profiles 

In our optical lattice, the spatial distribution of atoms 
is inhomogeneous because of the confining harmonic po- 
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FIG. 3: Metallic and insulating phases in the one-dimensional 
homogeneous superlattice system without harmonic confine- 
ment for Vsh/t — 6.0 and U/t — 4.0: (a) density of atoms, 
n, as a function of the chemical potential fi/t. The plateaus, 
which are formed clearly at n=0.5, 1.0, and 1.5, imply the 
existence of well-defined insulating phases, (b-l) and (b-2) 
are schematic diagrams of the insulating phases for n = 0.5 
and 1.5. 



FIG. 2: (a) Density profiles for a trap with N s = 120, N f = 
60, Vo/t = 0.01, U/t = 4.0 and V SL /t = 6.0. The filled circles 
denote the local density m and the solid lines are a guide for 
the eye. (b) shows the density profile, n" c , averaged within 
the unit cell. There are plateaus at n" c = 0.5, 1.0, and 1.5. 



tential. We thus consider the profile of the local density, 
i.e. the spatial distribution of rii, to be 



^(eta-Civ) 



(3) 



In Fig. Eta), we show the density profile for Nf = 60 
fermions trapped in a lattice with the linear size N s = 
120, where the other parameters are Vsh/t = 6.0, U/t = 
4.0 and V c /t — 0.01. It is seen in Fig. Et a ) that the am- 
plitude of the local density alternates between even and 
odd sites reflecting the 2-site periodic potential. When 
there is no superlattice structure, namely V"sl=0, it is 
known that a Mott insulating region appears character- 
ized by plateau formation at rii = 1.0 |2ll |22| . In the 
present case, however, the density profile reveals more 
complicated sawtoothed structures, from which we can 
not immediately distinguish metallic and insulating re- 
gions. To characterize the insulating region clearly, we 
therefore introduce the average number of atoms in the 
unit cell as, 



1 



-(rii + n l+1 ). 



(4) 



In Fig.Etb), we show the n" c value obtained by averaging 
the data of (a) in the unit cell. Plateaus (or shoulders) 
appear at n" c =0.5, 1.0, and 1.5, which characterize three 
distinct insulating regions. 

To understand the above characteristic behavior in 
nf c , it is instructive to observe the bulk properties of the 
system without harmonic confinement, which are shown 
in Fig. 3. Figure Et a ) shows the infinite DMRG results 
for filling factor n as a function of chemical potential 



\i. We can see the formation of well-defined plateaus at 
ri =0.5, 1.0, and 1.5, which clearly feature the insulating 
phases separated by the metallic phases. The insulating 
phases with n =0.5 and 1.5 are caused by the superlat- 
tice potential with 2-site periodicity. The corresponding 
spatial distribution of atoms is schematically illustrated 
in Fig. EHb). When n = 0.5, as shown in Fig. |3tb-l), the 
system is in a variant of the Mott insulating phase where 
each site with lower (higher) potential Vj accommodates 
one (no) atom. On the other hand, when n — 1.5, as 
shown in Fig. |3tb-2), the system is another Mott type 
insulator where each site with higher potential accom- 
modates one atom, while that with lower potential is al- 
most fully occupied by two atoms. The insulating phase 
with n = 1.0 changes its nature depending on the inter- 
action strength U in a slightly complicated way, from the 
band insulator (small U) to the Mott insulator (large U). 
Actually this problem has been extensively discussed so 
far, which is briefly summarized here. Early studies by 
the numerical 36] and weak-coupling bosonization meth- 
ods [33 suggested that a single phase transition occurs 
from the band insulator to the Mott insulator when U is 
varied at half filling. On the other hand, Fabrizio et al. 
predicted a new intermediate insulating phase "so-called 
bond-charge density wave (BCDW) phase" between them 
|38|. Since their proposal, extensive investigations have 
been done to confirm it. Quantum Monte Carlo simu- 
lations [13 supported the existence of the intermediate 
phase. By applying the method of topological transitions 
in spin and charge Berry phases, Torio et al. presented a 
ground-state phase diagram, which is consistent with the 
scenario of Fabrizio et al. 40] The DMRG calculations 
|4ll l42l I43L |44| have been performed by several groups. 
For example, the calculation of various structure factors 
supported the existence of intermediate phase 0] though 
Kampf et al. [SJ were not able to resolve two transi- 
tions clearly. Recent calculations based on the DMRG 
[45j and the level-spectroscopy method |46J with high ac- 
curacy have confirmed the existence of the intermediate 
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FIG. 4: (a) Profiles of the average density, n" c , along the 
trap for different fillings. The flat terraces are the insulating 
regions. The parameters are N 3 = 120, U/t = 4.0, Vsh/t = 
6.0, and V c /t = 0.005. In (a), the number of atoms is changed 
from 20 to 104, while in (b) the number is fixed as Nf =20, 
40, 60, and 90. 



BCDW phase. 

The above results for the homogeneous system show 
that the plateau regions at nf c —0.5, 1.0, and 1.5 in the 
system with confining potential have the above insulating 
properties, which are sandwiched by metallic regions. 

Let us now observe how the density profile of the sys- 
tem changes when the atom number Nf is changed sys- 
tematically. In Fig. Ufa), we show the average density 
of atoms, n" c , when Nf is changed continuously. It is 
seen that the system enters several distinct ground states 
as Nf increases. In Fig. Hfb), for reference, we show 
cross-sectional views of (a) by selecting several typical 
Nf values. For Nf — 20, there is an insulating region 
with nf c = 0.5 in the middle of the system. Just bed- 
side this region, there are narrow metallic region win- 
dows with nf c < 0.5, which smoothly continue to the 
vacant-atom region nf c ~ 0. The system is therefore in 
a phase where metals and insulators coexist, as is the 
case for an ordinary optical lattice without superlattice 



that the insulating regions with n" c = 0.5 are now sand- 
wiched by the two metallic regions with n" c < 0.5 and 
0.5 < n" c < 1.0. For Nf = 60, another insulating phase 
with n" c = 1.0 emerges in the middle region, and a fur- 
ther increase in the atom number, as seen when Nf = 90, 
results in an insulating region where n" c = 1.5. 

A brief comment on the finite-size effect is in order 
here. Since we are dealing with the rather small sys- 
tem with 60 fermions, the resulting metallic regions are 
very narrow, and the density profile changes very sharply 
there. This is due to the finite-size effect of our small 
system. If we consider a system with larger number of 
fermions, the metallic regions become wider and the den- 
sity profile changes more smoothly. 



B. Density fluctuations 

In order to characterize the above insulating regions, 
we further investigate local density fluctuations in our su- 
perlattice system. To this end, we introduce the variance 
of local density fluctuations, A,, at each site as 



A; 



(5) 



potentials |2lL |22|. For Nf = 40, a metallic region with 
0.5 < rC/ c < 1.0 appears in the middle of the system, so 



In ordinary one-dimensional optical lattices, it is known 
that Aj is suppressed and forms a plateau in the insu- 
lating region, while it has larger values in the metallic 
region |2lJ,|22j. This quantity, therefore, can describe the 
way in which insulator-metal crossover occurs when the 
on-site repulsion increases. 

Figure |SJ shows the local density n,;, the average density 
n" c , the variance Aj and the double-occupation proba- 
bility (ni^riii) for three different choices of U/t with the 
other parameters unchanged. 

It is convenient to discuss characteristic properties by 
specifying each region in terms of the average local den- 
sity n" c . Let us first observe the insulating region char- 
acterized by nf c = 1.0, which is located around the site 
i = 45 (or equivalently i — 75) in all three cases in Fig. El 
For U/t = 0.0, we can see the sawtoothed oscillation both 
in rii and (n^n^} since the system is in the band insula- 
tor where almost doubly-occupied sites and vacant sites 
are realized alternately. Note that even in this case, the 
variance Aj forms a plateau without such alternations, 
implying that density fluctuations occur locally only via 
a process between two adjacent sites. As the interac- 
tion strength U increases, m ((fJiT^iJ,)) continuously ap- 
proaches m — > 1 ((n^riii) 0). On the other hand, Aj 
is once enhanced (U/t = 3.0, Fig. E[b)), an d then sup- 
pressed (U/t = 9.0, Fig. Etc)). Therefore, we naturally 
expect that there may be another insulating state, which 
is characterized by an enhanced Aj, between the band 
insulating state (small U) and the Mott insulating state 
(large U). As mentioned above, in a homogeneous sys- 
tem without confinement, at half filling (corresponding 
to n l / c — 1.0), the BCDW intermediate phase appears 
between the band-insula ting phase and the Mott insu- 
lating phase lllii EI El El H El HI. In the 
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FIG. 5: Plots of the local density m, the average density n" c , the variance A, and the double-occupation probability (niftlii) 
for (a) U/t = 0.0, (b) 3.0, and (c) 9.0. The other parameters are fixed as N f = 60, iV s = 120, V c /t = 0.01, and Vsh/t = 5.0. 



BCDW state, there is rather large hopping within a cou- 
pled dimer, which gives rise to enhanced local density 
fluctuations. Therefore, we can distinguish the BCDW 
state from the other insulating states in terms of the vari- 
ance Aj. Summarizing the above results, we can say that 
the region with n" c = 1.0 is either in the band insulating 
state, the BCDW state or the Mott insulating state de- 
pending on the strength of the on-site repulsion. These 
states are smoothly connected to each other with the in- 
crease in U (crossover behavior) in contrast to the phase 
transitions realized in the homogeneous case. 

We next focus on the region around the site i = 35, 
where the average density n" c = 1.0 at U/t — 0.0 is 
changed to < c = 0.5 at U/t = 3.0 and 9.0. In this 
case, the local density rij is reduced from 2 to 1 at each 
odd site with the increase in U, and Aj and (jijfJijj.) 
decrease monotonically to zero (even sites are always 
empty). Namely, the system is smoothly driven from 
the band insulator to the Mott insulator. Finally, we 
look at the region around i = 60, where n" c is slightly 
larger than 1.5 at U/t = 0.0 without plateau structures, 
which is characteristic of the metallic region. In fact, 
Aj exhibits strong even-odd dependence on the site in- 



dex. At U/t — 3.0, we see the plateau formation at 
n" c = 1.5, where rij ~ 1.0 and {n^riii) ~ for even 
sites. Therefore, the region enters the insulating region 
with Mott (odd sites) and band-insulating (even sites) 
characteristics. A further increase in U again drives the 
Mott insulating state to the metallic state characterized 
by 1.0 < < c < 1.5. 

We have repeated similar calculations for various 
choices of the model parameters. Among them we have 
presented only the three selected cases above, since they 
capture most of the essential properties realized in our 
optical superlattice with 2-site periodicity. 



C. Spin correlations 

Finally, we discuss local and spatially-extended spin 
correlations in the system. In Fig. Ela)-(c), we show 
the computed results for the variance of local spin fluc- 
tuations ASi = ((Sf ) 2 ) together with the density pro- 
file rii and the average profile n" c for the system with 
V SL /t = 5.0 and U/t = 5.0. Also shown in Fig. 
(c) is the spin correlation function (SfS*) between the 
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FIG. 6: Plots of (a) the variance of local spin fluctuations 
ASi together with (b) rii and (c) n" c . The data are obtained 
for the system with Vsh/t = 5.0, U/t = 5.0. The other pa- 
rameters are fixed as N s = 120, N f = 60, V c /t = 0.01. 



FIG. 7: Plots of the spin correlation function between the ith 
and jth sites, (S*S$}, where the site j is fixed as (a)j = 35, 
(b)45 and (c)55. The data are obtained for the system with 
Vsh/t — 5.0, U/t = 5.0. The other parameters are fixed as 
N s = 120, N f = 60, V c /t = 0.01. 



ith and jth sites with the latter fixed as three typical 
positions, j = 35, 45, and 55. 

As clearly seen from Fig. \G\ the local moment is well 
developed (ASi > 0-2) when a given site is located in the 
Mott insulating region, while it is suppressed (ASi ~ 0) 
in the band-insulating region as well as in the empty 
region. These characteristic properties are in accordance 
with what we naively expect from the homogeneous case 
without harmonic confinement. Let us now look briefly 
at the spin correlation function. 

In Fig. 0a) , the i-dependence of the spin correlation 
(S*Sj) is shown where the site j — 35 is fixed in the 
insulating region with n" c = 0.5. It is seen that the spin 
correlation exhibits an antiferromagnetic nature within 
the same region, and quickly decays when the site i en- 
ters the metallic region. Interestingly, the spin correla- 
tion increases slightly again around i — 85 belonging to 
another insulating region with nf c = 0.5. This unusual 
behavior is discussed below in detail by using a more sig- 
nificant example. On the other hand, for j = 45 in the 
insulating region with n" c = 1.0, we cannot see the long- 
range spin correlation as in the j = 35 case, since this 
Mott insulating region is not well stabilized as regards 
the choice of parameters. For j = 55 in Fig. 0c), the 
spin correlation shows typical behavior inherent in the 
metallic state: a small spin correlation amplitude, which 
immediately decreases with distance. 

In Figs. 0D and ^\ we show the results obtained when 
we chose a slightly larger interaction strength U. In this 
case also, the spin correlation shows similar behavior to 
that in Fig. There are two points to be mentioned 
here. Since in this case the local moment is well devel- 
oped in the Mott insulating region, as seen from Fig. 
|5Ja), the spin correlation shows more significant antifer- 



romagnetic correlations. Even for site j belonging to the 
metallic region, we can observe slightly enhanced antifer- 
romagnetic correlations (Fig- Etc)). Another important 
point to be noticed is that the spin correlations in Fig. 
EJa) are strongly enhanced around i — 85, although they 
are suppressed in the region of i < 85 (c.f. Fig. 0a)). 
This characteristic behavior is a new feature that has not 
been observed in homogeneous systems without confine- 
ment. This is closely related to the formation of several 
distinct regions. In the Mott insulating region, a finite 
number of local spins are well developed, which are sensi- 
tive to effective magnetic fields. Such spins belonging to 
the different Mott regions can have enhanced correlations 
even if they are separated by the metallic phases. This 
is why we have encountered enhanced spin correlations 
between j = 35 and i — 85 in Fig. EUa). It should be 
noted that this effect is particularly significant if there 
is an odd number of spins in the Mott region. For ex- 
ample, in Figs. 0a) and^a), 3 sites around j = 35 and 
% = 85 belong to the well-defined Mott regions, which 
indeed exhibit enhanced spin correlations. 

Therefore, the above phenomenon is related to the 
finite-size effect of the Mott region, so that such behavior 
should be somewhat obscured when the Mott regions be- 
come large. It is thus interesting to know whether such 
behavior can really be observed by future experiments in 
the optical lattices. 



IV. SUMMARY 

We have studied the characteristic properties of 
fermions trapped in a one-dimensional optical superlat- 
tice with 2-site periodicity. We found that plateau re- 




FIG. 8: Similar plots to those in FigEJfor the other choice of FIG. 9: Similar plots to those in Fig|7]for the other choice of 
V SL /t = 3.0 and U/t = 7.0. Vsh/t = 3.0 and U/t = 7.0. 



gions appear in the profiles of the average number of 
atoms in the unit cell n" c =0.5, 1.0, and 1.5. In contrast 
to the ordinary Hubbard model in the optical lattice, new 
insulating regions appear with n" c =0.5 and 1.5 caused 
by the superlattice potentials. Furthermore, it has been 
shown by the variance and double-occupation probabil- 
ity results that three different insulating regions of band-, 
BCDW- and Mott-type emerge in the insulating region 
with n" c = 1.0, which are smoothly connected to each 
other via a crossover as the strength of the on-site repul- 
sion U is altered. 

We have also found that it is possible to enhance spin 
correlations between different Mott regions, even if they 
are once suppressed in the intermediate metallic region. 
This interesting phenomenon is due to the finite-size ef- 
fect of the insulating regions, which are naturally formed 
by harmonic confinement. Specifically, it occurs more sig- 
nificantly when an effective number of sites in the Mott 



insulating region is odd. 

Since systematic experimental studies on the one- 
dimensional optical superlattice may be soon within our 
reach, we hope that the characteristic properties dis- 
cussed for cold fermions in this paper will be observed 
in the near future. 
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